HRI-P-07-11-002 
CU-PHYSICS-12/2007 



X 



Unraveling neutrino parameters with a magical beta-beam 

experiment at INO 



Sanjib Kumar Agarwalla*'1^'", Sandhya Choubey* '', Amitava Raychaudhuri 

* Harish- Chandra Research Institute, 
Chhatnag Road, Jhunsi, Allahabad 211019, India 

'^Department of Physics, University of Calcutta, 
92 Acharya Prafulla Chandra Road, Kolkata 700009, India 



*,t,c 



ABSTRACT 

We expound in detail the physics reach of an experimental set-up in which the proposed large 
magnetized iron detector at the India-based Neutrino Observatory (INO) would serve as the far 
detector for a so-called beta-beam. If this pure z/g and/or z/g beam is shot from some source 
location like CERN such that the source-detector distance L ~ 7500 km, the impact of the CP 
phase 6c p on the oscillation probability and associated parameter correlation and degeneracies are 
almost negligible. This "magical" beta-beam experiment would have unprecedented sensitivity to 
the neutrino mass hierarchy and 613, two of the missing ingredients needed for our understanding 
of the neutrino sector. With Lorentz boost 7 = 650 and irrespective of the true value oi 6cp, the 
neutrino mass hierarchy could be determined at 3a C.L. if sin^ 2^13 (true) > 5.6 x 10~^ and we can 
expect an unambiguous signal for 613 at 3a C.L. if sin^ 26'i3(true) > 5.1 x 10"'^ independent of the 
true neutrino mass hierarchy. 
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1 Introduction 



Neutrino physics has entered the precision era, with the thrust now shifting to detailed under- 
standing of the structure of the neutrino mass matrix, accurate reconstruction of which would 
shed light on the underlying new physics that gives rise to neutrino mass and mixing. The full 
mass matrix is given in terms of nine parameters, the three neutrino masses, the three mixing 
angles and the three CP violating phases. Neutrino oscillation experiments are sensitive to only 
two mass squared differences, all the three mixing angles and the so-called Dirac CP violating 
phase. The remaining parameters, comprising of the absohitc neutrino mass scale and the two 
so-called Majorana phases, have to be determined elsewhere. We already have very good knowl- 
edge on the two mass squared differences and two of the three mixing angles. Results from solar 
neutrino experiments [1], which have been collecting data for more than four decades have now 
culminated in choosing the Large Mixing Angle (LMA) solution. The latest addition to this huge 
repertoire of experimental data is the result from the on-going Borexino experiment [2], and this 
result is consistent with the LMA solution. This conclusion from solar neutrino experiments has 
been corroborated independently by the KamLAND reactor antineutrino experiment [3, 4], and a 
combined analysis of the solar and KamLAND data gives as best-fit^ Am|i = 7.6 x 10~^ eV^ and 
sin^ 9i2 = 0.32 [4, 5]. The other mass squared difference Am|^ and mixing angle 623 are now pretty 
well determined by the zenith angle dependent atmospheric z/^ data in SuperKamiokande [6] and 
the long baseline experiments K2K [7] and MINOS [8]. The combined data from the atmospheric 
and long baseline experiments have pinned down lAmg^^l = 2.4 x 10~^ eV^ and sin^2^23 — 1- 

Despite these spectacular achievements, a lot of information is still required to complete our 
understanding of the neutrino sector. While the solar neutrino data have confirmed that Am2i > 
at a C.L. of more than 6a, we still do not know what is the sign of Am|]^. Knowing the ordering of 
the neutrino masses is of prime importance, because it dictates the structure of the neutrino mass 
matrix, and hence could give vital clues towards the underlying theory of neutrino masses and 
mixing. Knowing the sgnlAmli) could have other far-reaching phenomenological consequences. 
For instance, if it turns out the Am^^ < and yet neutrino-less double beta decay is not observed 
even in the very far future experiments, that would be a strong hint that the neutrinos are not 
Majorana particles (sec for e.g. [9] and references therein). Also, our knowledge on the third mixing 
angle 6*13 is restricted to an upper bound of sin^ ^^13 < 0.04 from the global analysis of all solar, 
atmospheric, long baseline and reactor data, including the CHOOZ [10] results in particular. Non- 
zero 9i3 brings in the possibility of large Earth matter effects [11, 12, 13] for GeV energy accelerator 
neutrinos travelling over long distances. Effect of matter on neutrino oscillations depends on the 
sgn{Amj^i) and is opposite for neutrinos and antineutrinos. For a given sgn{Amli) it enhances the 
oscillation probability in one of the channels and suppresses it in the other. Therefore, comparing 
the neutrino signal against the antineutrino signal in very long baseline experiments gives us a 
powerful tool to determine sgn{Aml{). A non-zero value of this mixing angle would also open up 
the possibility of detecting CP violation in the neutrino sector. 

Tremendous effort is underway to determine 613, the CP phase dcp and sgn{Am^^) using long 
baseline experiments [14, 15, 16, 17]. Future programs involving accelerator based neutrino beams 

^We use a convention where Amf^ = m? — m^. 
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include among others the T2K [14], NOz/A [15], Superbeams, beta-beams and Neutrino Factories 
[17]. In earher papers [18, 19] we have proposed and expounded the possibihty of measuring 
to a very high degree of accuracy the mixing angle ^13 and sgnlAmli) aka, the neutrino mass 
hierarchy^, in an experimental set-up where a pure and intense Ue and/or z/g beam is shot from 
CERN to the India-based Neutrino Observatory (INO) [20]. This pure and intense source of i^e 
and/or flux could be the so-called beta-beam [21], which is created when fully ionized and highly 
accelerated radioactive ions beta decay in the straight sections of a ring, where they are circulated 
and stored, after being produced, collected, bunched and accelerated. A large magnetized iron 
calorimeter (ICAL) is expected to comc-up soon at the INO facility in India. Since the energy 
threshold of this detector would be at least 1 GcV, the beta-beam should necessarily be a multi- 
GeV beam. While the most widely discussed source ions for beta-beams, ^^Ne and ^He, need very 
large values of the Lorentz boost to reach multi-GeV energies [18], alternative ions with larger 
end-point energy, such as^B and ^Li, can be used with reasonable acceleration. While the most 
discussed design for the beta-beam set-up [22] needs modest Lorentz boosts, it suffers from the 
effect of the so-called "parameter degeneracies" [23, 24, 25] giving rise to eight-fold degenerate 
"clone solutions" [26]. A very big advantage that the CERN-INO beta-beam experiment would 
have is that at the baseline of 7152 km, the Sep dependent terms (almost) drop out from the 
expression for the — > oscillation probability. As a result, this experiment is free of two of the 
three degeneracies [26, 27, 28]. In addition, large energies and the large distance involved allows 
the neutrino to pick up near-resonant matter effects, enhancing the oscillation probability which 
can thereby compensate to a great extent the reduction of the neutrino flux due to the 1/L^ factor 
[19]. This makes the CERN-INO beta-beam set-up almost magical for determining sgn{Amli) 
and ^13. 

In [19] we studied the physics potential of this experimental set-up when we run the beta-beam 
in only one polarity for five or ten years. That is, wc probed the sensitivity of the experiment to 
sgn{Am^i) and using only the Ue {ve) beam for five years running with 1.1 x 10^^ (2.9 x 10^^) 
useful ion decays per year. In this paper we extend our analysis by including data from both the 
neutrino and antineutrino run of this beta-beam experimental set-up. We demonstrate how adding 
data from both polarities serves to strengthen precisely those regions where the individual ones are 
less powerful and thus significantly enhances the mass hierarchy sensitivity of the experiment. As 
a further refinement, we analyse the full spectral data expected in the CERN-INO beta-beam set- 
up. We present and compare results from the rates-only analysis against results from the analysis 
where we incorporate the full spectrum. For ^13 measurement, we consider two scenarios: (i) when 
there are no (or v^) events in the detector and (ii) when we see a signal in the detector. In the 
former case we present the 3o" upper bound on sin^ 2^^i3 expected from the null results. We call this 
the "sin^ 26*13 sensitivity reach" of the experiment. In the latter case, we first study the range of 
"true" values^ of sin^ 26'i3(true) for which the experiment would be able claim to have seen a signal 
at the 3(7 C.L. This is termed as the "sin^ 2^13 (true) discovery reach" of the experiment. Finally, 
for a non-zero signal at the detector, we study how precisely sin^ 2^13 can be determined with 5 

^What we usually refer to as the neutrino mass hierarchy is really the neutrino mass ordering. Therefore, our 
discussion and results are equally relevant for a quasi-degenerate neutrino mass spectrum as they are for hierarchical 
and inverted hierarchical spectra. 

^Throughout this paper we denote true value of the parameters by putting "(true)" in front of them. 
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Figure 1: Left panel shows the total flux in yr~^m~'^ expected at INO, as a function of the 
Lorentz factor 7. The solid (dashed) line corresponds to i^Li) and we have assumed 1.1 x 10^*^ 
(2.9 X 10^^) useful ions decays per year. Right panel shows the energy at which the flux peaks, as 
a function of 7. 



years of combined neutrino and antineutrino run. In our earlier paper [19], we had given all results 
assuming that 5cp (true) = 0. The hierarchy and ^13 measurement sensitivities however depend 
on the value of 5cp (true). In this paper we find the physics reach of the CERN-INO beta-beam 
set-up for all possible values of 5c'p(true) = and show the best and worst possible physics reach. 
We also study the impact of changing the number of useful ion decays in the straight sections of 
the beta-beam storage ring and compare the dependence of the sensitivity on the number of ion 
decays and the Lorentz boost 7. 

For all results presented in this paper, we use the full PREM Earth matter density profile 
[30] for simulating the prospective data. When we fit this simulated data, we allow for a 5% 
uncertainty in the PREM profile and take it into account by inserting a prior and marginalizing 
over the density normalization. We also study the impact of changing the Earth matter density by 
±5% in the data itself. We also study the impact of changing the energy threshold of the detector 
and the background rejection factor. In our analysis here we also allow the parameters Am^i and 
sin^ 612 to vary arbitrarily in the fit. 

We begin by providing a brief overview of the proposed experimental set-up, the expected 
event rate, and the importance of the magic baseline in section 2. In section 3, we give our results 
on the sensitivity of the CERN-INO beta-beam set-up to the neutrino mass hierarchy. In section 

4, we discuss the potential of measuring/constraining ^13. We end with our conclusions in section 

5. The details of our numerical code and analysis procedure are relegated to an Appendix. 
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Figure 2: Total number of expected events in five years as a function of the baseline L for the 
source with 7 = 500 and for two values of sin^ 26'i3 and assuming that the normal hierarchy is 
true. The hatched areas show the expected uncertainty due to the CP phase. 

2 The CERN-INO beta-beam Experimental Set-up 

Very pure and intense Pe and/or z/g beams can be produced by the decay of highly accelerated 
radioactive beta unstable ions, circulating in a storage ring. This is what is called a "beta-beam" 
and was first proposed by Piero Zucchelli [21]. Since the flux spectrum is determined entirely by 
the end-point energy of the parent ions and the Lorentz boost provided by the acceleration, it is 
almost free of systematic uncertainties. The flux normalization is determined by the number of 
useful ion decays in the straight section of the beam. The selection of the beta unstable parent 
ion is determined by a variety of factors essential for efficiently producing, bunching, accelerating 
and storing these ions in the storage ring. Among the candidate ions discussed in the literature, 
the ones which have received most attention so far are the ^^A^e and ^He ions for producing the Ve 
and z/g beams respectively [22, 31, 32, 33, 34, 35]. Both these ions have a similar end-point energy 
which is about 4 MeV in the parent ion's rest frame. Two other candidates, and ^Li, as source 
ions for z/g and z/g beams respectively, have been recently shown to be viable [36, 37, 38, 39, 40]. 
The advantage that these ions have over ^^Ne and ^He is their larger end-point energies, which 
is higher by a factor of more than 3. Therefore, for the same Lorentz boost factor, we expect 
the resultant ^B and ^Li beams to be about 3 times higher in energy compared to the ^^A^e and 
^He beams. We refer the reader to Table 1 of our earlier paper [19] for the full details about 
characteristics of the four beta-beam candidate ions. 

These large beta decay end-point energy ions are particularly important for the CERN-INO 



5 




Figure 3: Total number of events as a function of 7 for the source, for different values of L 
are shown in the four panels. The black hatched area shows the uncertainty range due to the CP 
phase when NH is true, while the area between the maroon dashed lines shows the corresponding 
uncertainty when IH is true. For all cases we assume sin^ 2^13 = 0.05. 



beta-beam experimental set-up discussed in [18, 19], where the beta-beam from CERN is shot 
to a magnetized iron detector at the India-based neutrino observatory (INO). The INO will be 
located in southern India, close to the city of Bangalore. The CERN to INO distance corresponds 
to 7152 km, which is tantalizingly close to the "magic baseline" [27, 28] (see also [29]). Being 
free of the problem of parameter degeneracies [23, 24, 25, 26], the magic baseline is known to be 
particularly useful for measuring the mixing angle 6'i3 and sgnlArn^i). The concept of the magic 
baseline can be very easily understood by looking at the approximate expression of Pe/i? where 
the conversion probability is expanded in the small parameters 6*13 and a = Am^i/Amg^, keeping 
only terms up to second order in these small parameters. This expression for the "golden channel" 
[41] probability is given as [41, 42] 

o n sin2[(l-i)A] 
Pen ^ sin^ ^23 sin^ 2^13 ^ 

sin(iA) sin[(l - i)A] 



-f a sin 26'i3 sin 2^12 sin 26'23 sin 6cp sin(A) 



A {1-A) 
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Figure 4: Total number of events as a function of 7 for the source, for different values of L are 
shown in the four panels. The black hatched area shows the uncertainty range in the events due to 
CP phase when sin^26'i3 = 0.05, while the red hatched area shows the corresponding uncertainty 
when sin^2^i3 = 0.01. For all cases we assume NH to be true. 



. . Nsin(iA) sin[(l - i)Al 
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l\m%L A 



where A = ■^2\/2GpNf.E is the matter potential, N^, being the electron number density inside the 
earth and Gp the Fermi constant, the + sign refers to neutrinos while the — to antineutrinos. In 
Eq. (1) the second term has the CP violating part. The CP phase 5c p is positive for neutrinos and 
negative for antineutrinos and therefore the sin^cp term changes sign. The third term, though 
5c p dependent, is CP conserving, while the fourth term is independent of ^13 as well as 5c p- If 
there exists a baseline L for which the condition 

sin(iA) = (3) 
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Figure 5: Total number of events as a function of 7 for the (solid lines) and the ^Li (dashed 
lines) sources. Results for both normal and inverted hierarchies are shown. 



holds, then the second, third and fourth terms in Eq. (1) drop out and one is left with just the 
first term. In particular, we see that all 6cp dependent terms go away at this magic baseline. 
Therefore, this experimental set-up is free of two of the three parameter degeneracies, providing 
us with a firm bedrock for determining ^13 and the mass hierarchy [27, 19]. It turns out that for 
the PREM Earth matter density profile, the condition given by Eq. (3) is satisfied for L ~ 7500 
km. In [19] we also stressed the point that for this very long baseline, neutrinos would also pick 
up large, and possibly near-resonant, matter effects. Largest enhancement of oscillations due to 
matter effects of course comes about when the product of the mixing angle term in matter and 
the (Amli)'^^ driven oscillatory term is the largest [43, 44]. For the CERN-INO baseline of 7152 
km, the probability for largest conversion is expected for E 6 GeV, for sin^ 2^13 = 0.01 and 
Amii = 2.5 X 10-3 eV^. 

The ICAL detector, which will be a 50 kton magnetized iron calorimeter, will be built at INO 
[20]. There is a possibility that the detector mass might be increased to 100 kton at a later stage. 
The approved design of the detector comprises of 6 cm iron slabs interleaved with 2 cm thick glass 
Resistive Plate Chambers (RPC), which would serve as the active detector elements for ICAL. 
The iron will be magnetized by an external field of about 1 Tesla, giving the detector charge 
identification capability. The detection efficiency of ICAL after cuts is expected to be about 80% 
and energy threshold would be about 1 GeV. In what follows, we will use an energy threshold of 
1.5 GeV for our main results. However, we will also show the impact of changing the threshold. 
The energy resolution of the detector is expected to be reasonable and we assume that the neutrino 
energy will be reconstructed with an uncertainty parameterized by a Gaussian energy resolution 
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function with a HWHM aE = 0.15E, where E is the energy of the neutrino. We will present and 
compare the sensitivity of this experimental set-up with and without the full spectral analysis. 
Details of our numerical approach can be found in [19] and in the Appendix. 

In the left panel of Fig. 1 we show the energy integrated total number of neutrinos in units of 
yr~^m~'^ arriving at INO, as a function of the Lorentz boost 7. The solid (dashed) line corresponds 
to (^Li) and we have assumed 1.1 x 10^^ (2.9 x 10^^) useful ion decays per year*^. Throughout 
this paper we have assumed that the detector is aligned along the axis of the beam. The figure 
shows that the energy integrated flux arriving at the detector increases almost quadratically with 
7. Note that with the same accelerator, the Lorentz boost acquired by is 1.67 times larger than 
that by ^Li, determined by the charge to mass ratios of the ions. The right panel of Fig. 1 depicts 
the energy at which the flux peaks, as a function of the Lorentz boost 7. It turns out that this peak 
energy is roughly half the maximum energy of the beam, which is given as E'max — 2 (£'0 — me) 7. 

In Fig. 2 we show the number of events expected in five years as a function of the baseline L, 
if wc run the experiment in the neutrino mode with 7 = 500. A similar figure is expected for the 
antincutrino beam. The upper black hatched area shows the events for sin^ 2^13 = 0.05 and the 
lower red hatched area corresponds to sin^ 2^13 = 0.01. For each baseline L, the range covered 
by the hatched area shows the uncertainty in the expected value of the number of events due to 
the completely unknown dcp, which could take any value from to 27r. The baseline L where 
the width of this band reduces to (almost) zero is the magic baseline. We see from the figure 
that the magic baseline appears at about L ~ 7500 km. Note that while for sin^ 2^^i3 = 0.01 the 
magic baseline is very clearly defined with the CP dependence going completely to zero, for the 
higher value of sin^ 2^13 of 0.05, the "magic" is not complete. The reason for this anomaly can be 
traced to the fact that Eq. (1) was derived for only very small values of ^13. For larger values of 
this angle, higher order terms become important. These terms might depend on 5cp and remain 
non-zero even at the magic baseline. 

In Figs. 3 and 4 we show the impact of the magic baseline on the mass hierarchy and 6*13 
sensitivity respectively. In each of the four panels of both the figures we show the expected 
events in five years as a function of 7. Each panel is for a certain fixed value of L, shown in 
the corresponding panel. In Fig. 3 the black hatched area shows the band for normal hierarchy 
(NH) while the open band delimited by the dashed red lines are for the inverted hierarchy (IH). 
As in Fig. 2 the band correspond to the uncertainty in the event rate due to the unknown 5cp- 
The effect of the uncertainty of 5c p almost vanishes for L — 7500 km which is very close to the 
magic basehne. We can see that for the smaller baseline L — 1000 km, NH and IH predictions are 
largely overlapping, making it almost impossible for these experiments to give sensitivity to the 
mass hierarchy unless sin^ 2^^i3 turns out to be extremely large and 5cp favorable. The hierarchy 
sensitivity is expected to improve as we go to larger baselines and this is reflected from the two 
bands for NH and IH separating out. It turns out that because the matter effects are very large 
for the magic basehne and effect of CP uncertainty is zero, this basehne gives the best sensitivity 
to the mass hierarchy. For L larger than magic, matter effects are higher but the flux is lower, 
while for L lower than magic, flux is higher but the matter effects are lower. For both above and 

^Unless stated otherwise these are the reference luminosities in all the figures. Also, all figures correspond to a 
five year run. 
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below the magic baseline, the effect of 6cp is expected to further reduce the sensitivity. This is 
particularly true for the lower L baselines. Fig. 4 shows the bands for NH but with two different 
choices for sin^ 2^13. Here the effect of the magic basehne is seen even more clearly. 

In Fig. 5 we show as a function of 7, the number of events expected in five years in the 
CERN-INO beta-beam set-up. The solid (dashed) lines are for neutrino (antineutrino) events, 
with the thick line showing the event rate for NH (IH) while the thin line is for the IH (NH). We 
have assumed sin^ 2^13 = 0.01 and 5cp — 0. One point which is transparent from this figure and 
which will be very relevant in understanding the behavior of the CERN-INO beta-beam set-up 
is the following: For a given value of ^13 and for NH (IH), we expect a large number of events 
in the neutrino (antineutrino) channel and almost negligible events in the antineutrino (neutrino) 
channel. This means that for NH (IH) it will be the neutrino (antineutrino) channel which will 
be statistically more important. 

As discussed in detail in [19] we expect hardly any background events in the CERN-INO beta- 
beam experiment. The atmospheric neutrino flux falls steeply with energy and is expected to 
produce much fewer events for the energy range that we are interested in^. The fact that INO has 
charge identification capability further reduces the atmospheric background. The most important 
handle on the reduction of this background comes from the timing information of the ion bunches 
inside the storage ring. For 5T magnetic field and 7 = 650 for ions, the total length of the 
storage ring turns out to be 19564 m. We have checked that with eight bunches inside this ring at 
any given time, a bunch size of about 40 ns would give an atmospheric background to signal ratio 
of about 10~^, even for a very low sin^ 26'i3 of 10~^. For a smaller bunch span, this will go down 
even further. In addition, atmospheric neutrinos will be measured in INO during deadtime and 
this can also be used to subtract them out. Hence we do not include this negligible background 
here. In our numerical analysis we have assumed that the background events come only from 
the neutral current interactions of the beta-beam in the detector. Wc estimate it by assuming 
an energy independent background fraction of ~ 10"*^ [45]. We have noted that after five years 
of running of the CERN-INO beta-beam experiment with 7 = 650, we expect only about 0.1 
background events. Nevertheless we take this background into account in our numerical analysis. 
Details of our numerical analysis are given in the Appendix. 

3 Measurement of the Neutrino Mass Hierarchy 

In Fig. 6 we show the minimum value of sin^ 2^13 for which it would still be possible to reject the 
wrong hierarchy by this experimental set-up at the 3a C.L.^. The left panel shows the hierarchy 
sensitivity when the normal hierarchy is true, while the right panel corresponds to the inverted 
hierarchy being true. We show this as a function of 7. In both panels we show by the red solid 
lines the sensitivity when we add neutrino and antineutrino data with the same value of the 
Lorentz boost shown in the x-axis. The blue dashed lines on the other hand correspond to the 
sensitivity expected when the neutrino beam runs with the 7 shown in the x-axis while the 7 for 

^We will show in the next section that even a threshold energy of 4 GeV is easily admissible in our set-up, and 
above 4 GeV there are much fewer atmospheric events. 

^We have given the details of our numerical analysis method in the Appendix. 
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Figure 6: Minimum value of sin^ 2^13 (true) for which the wrong hierarchy can be ruled out at the 
3(7 C.L., as a function of 7. The left panel is for normal hierarchy as true, while the right panel is 
when inverted hierarchy is true. The red solid curves show the sensitivity when the 7 is chosen to 
be the same for both the neutrino and the antineutrino beams. The blue dashed lines show the 
corresponding sensitivity when the 7 for the antineutrinos is scaled down by a factor of 1.67 with 
respect to the 7 of the neutrino beam. 



the antineutrino beam is scaled down by a factor of 1.67. We assume five years of running of 
the beta-beam in both polarities and we have performed a full spectral analysis. We note that 
using the combined neutrino and antineutrino beam running at the same value of 7 significantly 
enhances the hierarchy sensitivity of the experiment and the wrong hierarchy could be ruled out at 
3(7 for sin^ 26'i3 > 6.8 x 10~^ (sin^ 2^13 > 4.0 x 10~^) for 7 = 500 if the normal (inverted) hierarchy 
is true. This should be compared to the results from our earlier paper where we demonstrated 
that by using the beam with a single polarity, the wrong hierarchy can be rejected at the 3a 
CL. for sin^ 26'i3 > 9.8 x 10"'^ (sin^ 2^13 > 8.2 x 10~^) when the normal (inverted) hierarchy is 
true. Partial improvement in the hierarchy sensitivity comes from the fact that we have used 
spectral information here, while in our earlier analysis we had considered just the total event rate 
measured in the detector. However, here the major improvement comes due to the addition of 
both the neutrino and antineutrino events, while in [19] we had considered events due to either 
the neutrino (for true normal hierarchy) or the antineutrino (for true inverted hierarchy) beam 
alone. Presence of both neutrino and antineutrino data simultaneously in the analysis restricts 
the fitted value of ^13 to be in a range very close to the assumed true value. For instance, for NH 
true, data corresponds to a large number of events for neutrinos and a small number of events 
for antineutrinos. When this is fitted with IH, we have a small number of events predicted for 
the neutrinos. In order to minimize the disparity between the data and prediction for neutrinos, 
the fit tends to drive 6^13 to its largest allowed value. However, larger values of 6^13 would give 
very large number of antineutrino events for IH and this would be in clear conflict with the data. 
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Figure 7: Plots showing the impact of various factors on the mass hierarchy sensitivity of the 
CERN-INO beta-beam experiment. The top left panel shows the impact of changing the detector 
threshold. The lower left panel shows the effect of changing the background rejection factor. The 
top right panel shows the difference in the sensitivity between the rate and spectral analysis. The 
lower right panel shows how the density profile would impact the hierarchy sensitivity. 
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Figure 8: Effect of 6cp{tTue) on tlie fiierarchy sensitivity. Tlie black daslied lines show the worst 
and best cases when we allow ^c-p (true) to take any value between and 27r. The red solid curve 
corresponds to the reference case where 6cp{tTue) = 0. The left panel shows the case for true 
normal hierarchy while the right panel is for true inverted hierarchy. 



Therefore, the net advantage of adding data from both neutrino and antineutrino channels is that 
one cannot artificially reduce the ^iny longer by tinkering with ^13 in the fit. As a result, the 
sensitivity of the experiment to mass hierarchy witnesses a substantial improvement. 

We note from the plots that the hierarchy sensitivity falls when we use the scaled 7 option for 
the antineutrino beam. This is particularly relevant when the true hierarchy is inverted and/or 
when 7 is low. Since scaling the 7 reduces it by a factor of 1.67, the statistics for the antineutrinos 
fall by nearly a factor of 1.67 for this case and this reflects in the reduced hierarchy sensitivity of 
the experiment. Its impact when true hierarchy is inverted is more because in that case, the data 
corresponds to larger events for the antineutrinos and very small events for the neutrinos. The 
antineutrino events are therefore the driving force and an increase in their statistical uncertainty 
due to the scaled down 7 accentuates the adverse effect on the hierarchy sensitivity. In the case 
of normal hierarchy, the events in the neutrino channel are the dominant factor and the role 
of the antineutrinos is only to prevent the 613 values in the fit to run to very large values, as 
discussed before. As long as the antineutrino events have enough statistical power to restrict ^^13 
to values close to the true value at which the data was generated, the hierarchy sensitivity remains 
reasonably good. Therefore, for the normal hierarchy only for very low values of 7 the hierarchy 
sensitivity gets seriously affected by the Lorentz boost scaling. 

In Fig. 7 we show how the hierarchy sensitivity depends on diverse input factors. As in Fig. 
6 we show the 3a limit for sin^ 2^13 as a function of 7 in all the four panels and we assume that 
normal hierarchy is true. The reference curve (red solid line) in all panels corresponds to the 
result obtained with a z/g and z/g beam with a spectral analysis. The upper left hand panel shows 
the effect of changing the threshold energy of the detector. The sensitivity of the experiment is 
seen to remain almost stable against the variation of the threshold energy upto 4 GeV. Only for a 
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Figure 9: The variation of the experimental sensitivity on the number of useful ion decays in the 
straight sections of the storage ring. Left panel shows sensitivity to the mass hierarchy assuming 
NH to be true. Right panel shows the sin^ 2^13 sensitivity reach. In both panels, the majenta solid 
vertical line corresponds to the reference value used in the rest of the analysis. 



threshold of 6 GeV and above the sensitivity falls, the lower 7 values getting more affected since 
they correspond to lower neutrino energies. In the lower left hand panel of the figure we show 
the effect of the chosen background fraction on the hierarchy sensitivity. The red solid line shows 
the sensitivity for our assumed background factor of 10~^ while the black dashed line shows the 
corresponding sensitivity when the background rejection is poorer and we have a higher residual 
background fraction of 10~^. The upper right hand panel shows how our sensitivity increases by 
taking into account the spectral information of the events. It also shows how much improvement 
we get by combining the antineutrino data with the neutrino data. The black dashed line shows 
how the sensitivity falls when we use the total event rates instead of the events spectrum. The 
blue dashed-dotted and green dashed-triple-dotted hues show the sensitivity expected from the 
neutrino data alone. The blue dashed-dotted line is for binned neutrino data while the green 
dashed-triple-dotted lines shows the sensitivity for the total event rate for neutrinos alone. We 
see that the effect of using the spectral information is only marginal when both neutrino and 
antineutrino are used, while the effect of combining the antineutrino data with the neutrino data 
on the sensitivity is huge. For the neutrino data alone, the sensitivity improves significantly when 
one uses the spectral information. In the lower right hand panel we show how the sensitivity of 
the experiment to hierarchy would get affected if we use a different profile for the Earth matter 
density instead of PREM. The red solid line is for earth density according to the PREM profile 
while the blue dotted and black dashed lines are when the matter density is 5% lower and 5% 
higher respectively than the density predicted by the PREM profile. When the density is higher 
(lower) the matter effects are higher (lower) and therefore the sensitivity improves (deteriorates). 

One crucial point that we have not stressed so far concerns the dependence of the detec- 
tor performance on the true value of 6c p- AH the earlier plots were presented assuming that 
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7\N 


Mass Hierarchy (3o") 


sin^ 29i^ sensitivity (3o") 


1.1 X W 


2.043 X 10^» 


1.1 X W 


2.043 X W 


350 
650 


1.3 X 10-3 
5.6 X 10-^ 


9.3 X 10"^ 
4.1 X 10-^ 


3.8 X 10-3 
1.1 X 10-3 


2.3 X 10-3 
7.3 X 10-^ 



Table 1: Comparison of the variation of the detector sensitivity to mass hierarchy (columns 2 and 
3) and sin^ 26'i3 sensitivity (columns 4 and 5) with 7 and N, the number of useful ion decays per 
year. 



Scp{tTne) = 0. At exactly the magic baseline, we expect the sensitivity of the experiment to be 
completely independent of 6cp- The CERN-INO distance of 7152 km is almost magical, but it 
is not the exact magic baseline. Therefore, we do expect some remnant impact of 6cp{tTue) on 
our results^. To show how our results get affected by 5cp{tTue), we show in Fig. 8 the hierar- 
chy sensitivity just as in Fig. 6, but here we show the full band corresponding to all values of 
Scp{tTue) from to 27r. As before, the left panel is for NH true while the the right panel is for 
IH true, and we have taken in the analysis the full spectral data for the neutrinos as well as the 
antineutrinos, with the same 7. The lower edge of this band shows the best possible scenario 
where the experiment is most sensitive, while the upper edge shows the worst possible sensitivity. 
The red solid lines in both panels show for comparison the hierarchy sensitivity corresponding 
to 5cp(true) = 0, which we had presented in Fig. 6. We note from the figure that the hier- 
archy sensitivity is nearly the best for 6cp{tTne) = when IH is true while if NH is true then 
it would give us almost the worst sensitivity. For NH (IH) as true the best possible sensitivity 
would be 8111^2913 > 3.96 x 10"^ (sin2 2^i3 > 2.96 x 10"^) for 7 = 650 to be compared with 
sin2 2^i3 > 5.51 x 10"^ {8111^2913 > 3.05 x 10"^) when Scp{true) = 0. Therefore, we conclude 
that if NH is true then it would not be unfair to expect an even better hierarchy sensitivity than 
what was reported in Fig. 6, while if IH is true then the best sensitivity will be returned for 
5cp(true) ~ 0. 

We have noted from Figs. 3 and 4 that the total number of events in the detector increases 
roughly linearly with 7, except for extremely long baselines. Increasing the number of ion decays 
per year will also bring about a simple linear increase in the statistics. It is therefore pertinent to 
make a fair comparison between the dependence of the mass hierarchy sensitivity to the Lorentz 
boost 7 and the number of useful ion decays in the ring®. In the left panel of Fig. 9 we show the 
effect of increasing the number of ion decays on the hierarchy sensitivity^. The plots exhibit the 

''Note that in all our results presented in this paper, we have fully marginalized over all the oscillation parameters 
in the fit, including 6c p- 

^Note that this is also equivalent to increasing the total exposure time of the experiment. Both number of ion 
decays per year and exposure appear as a normalization factor for the event rate and hence increasing the number 
of ion decays by a factor n keeping the exposure same is equivalent to increasing the exposure by a factor n keeping 
the number of ion decays per year fixed. 

^We assume that the number of useful ion decays for both and ^Li have been scaled by the same factor. In 
the figure along the a;-ax;is only the numbers are shown. 
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Figure 10: Left panel shows the 3a sensitivity hmit for sin^ 26'i3. Right panel shows the 3a discovery 
reach for sin^ 2^13 (true). The red solid lines in the left and right panels show the sensitivity reach 
and discovery potential respectively, when the 7 is assumed to be the same for both the neutrino 
and the antineutrino beams. The blue dashed lines show the corresponding limits when the 7 for 
the ^Li is scaled down by a factor of 1.67 with respect to the 7 of the neutrino beam, which is 
plotted in the x-axis. 

dependence of the sensitivity on the number of useful ion decays per year for an exposure of five 
years, for three different values of 7. We have assumed the same Lorentz boost for the neutrino 
and antineutrino beams. We present in Table 1 the relative increase in the hierarchy sensitivity 
when we increase the 7 by a factor of 1.86 and compare it against the increase in the sensitivity 
when the number of ion decays are increased by the same factor. We note that while the hierarchy 
sensitivity improves by a factor of 2.54 in going from 7 = 350 to 650 keeping the number of ion 
decays per year as 1.1 x 10^^, it increases 1.5-fold when we raise the number of ion decays per 
year from 1.1 x 10^^ to 2.04 x 10^^ keeping 7 = 350. However, we would like to stress that the 
improvement of the hierarchy sensitivity is not linear with either 7 or number of ion decays per 
year. The crucial thing is that the behavior of the sensitivity dependence on both 7 and number of 
ion decays per year is very similar. It increases very fast initially and then comparatively flattens 
out. 

4 Measurement of sin^ 2^13 

The CERN-INO beta-beam set-up is also expected to give very good sensitivity to the 613 mea- 
surement. In what follows, we will quantify our results in terms of three "performance indicators", 

1. sin^ 26'i3 sensitivity reach, 

2. sin^ 26'i3 discovery reach. 
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3. sin^ 29i3 precision. 

We give below a detailed description of our definition of these performance indicators. For all 
results in this section we take into account the full event spectrum and combine five yeats data 
from both the neutrino and antineutrino channels. 

4.1 sin^2^i3 Sensitivity Reach 

We define the sensitivity reach of the CERN-INO beta-beam experiment as the upper limit on 
sin^ 26*13 that can be put at the 3a C.L., in case no signal for ^^13 driven oscillations is observed 
and the data is consistent with the null hypothesis. We simulate this situation in our analysis by 
generating the data at sin^ 2^13 (true) = and fitting it with some non-zero value of sin^2^i3 by 
means of the technique. In our fit we marginalize over all the oscillation parameters including 
Sep and further marginalize over the mass hicrarchy^'^ and choose the value of sin^ 2^^i3 for which 
the fit yields = 9- The result is shown in the left panel of Fig. 10, as a function of 7. The red 
solid line shows the sin^ 2^13 sensitivity when 7 is assumed to be the same for both the neutrino 
and the antineutrino beams. The bhic dashed lines shows the corresponding 3a upper limit when 
7 for the ^Li is scaled down by a factor of 1.67 with respect to that for the neutrino beam. We 
can compare our results here with those obtained from our earlier analysis in [19] where we had 
considered data from either the neutrino or the antineutrino channel alone when the hierarchy was 
normal and inverted, respectively. We find that the sin^ 2^13 we obtain from the combined neutrino 
and antineutrino data does not exhibit any marked improvement compared to that obtained in 
[19]. This can be understood by the following reasoning. In generating the data we have assumed 
that sin^ 2^13 (true) = 0, which means that we have negligible events in both the neutrino as 
well as the antineutrino channels, irrespective of the mass hierarchy. When this data is fitted 
allowing for non-zero sin^2^i3, the neutrino (antineutrino) channel plays a dominating role when 
NH (IH) is assumed in the fit. Therefore, the sensitivity we obtain assuming NH (IH) in our fit is 
similar to what we had in [19] for the neutrino (antineutrino) channel alone. However, we reiterate 
that Fig. 10 shows the sin^ 2^^i3 sensitivity after marginalizing over hierarchy as well. In other 
words, the sensitivity shown in this figure corresponds to the statistically weaker channel. For the 
case where we use same 7 for and ^Li, the neutrino channel is weaker since the event rate is 
about 1.5 times less than antineutrino events with the same 7. On the other hand when we scale 
down the Lorentz boost for ^Li, the fiux in the antineutrino channel goes down significantly and 
hence it becomes the statistically weaker channel as can be seen from Fig. 5 and therefore the 
marginalized corresponds mainly to that from antineutrinos. Indeed one can check that the 
sin^ 2^13 sensitivity that we exhibit by the blue dashed line for the scaled 7 case is comparable to 
what we had obtained for the antineutrino channel with IH and the corresponding lower 7. 

The dependence of the sin^ 2^13 sensitivity on the number of useful radioactive ion decays per 
year in the straight section of the storage ring is shown in the right panel of Fig. 9. Here we 
have taken the same Lorentz boost for and ^Li and we have shown the results for three fixed 

^°Note that since sin^ 2^13 (true) = 0, the data is independent of the mass hierarchy. However, since we allow for 
non-zero sin^ 2^13 in the fit, the predicted event rates in our "theory" depend on the mass hierarchy. 
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values of 7. The relative increase in the sensitivity by increasing 7 and/or the number of useful 
ion decays per year by the same factor is quantified in the last two columns of Table 1. 

4.2 sin^ 2^13 Discovery Reach 

How good are our chances of observing a positive signal for oscillations and hence 6*13 in the CERN- 
INO beta-beam set-up? We answer this question in terms of the parameter indicator which we call 
the "discovery reach" of the experiment for sin^ 2^13. This is the minimum value of sin^ 26*13 (true) 
that would give an unambiguous signal in the detector at So". To find this we simulate the data at 
some non-zero value of sin^ 26^13 (true) and fit it by assuming that sin^ 26^13 = 0, allowing all other 
oscillations parameters to take any possible value in order to return back the smallest value for 
the x^. Note that since the fitted value of the mixing angle in this case always corresponds to 0, 
there is no need of any marginalizing over the hierarchy when fitting the data. However, since the 
data here is generated at a non-zero value of sin^ 26^13 (true), it depends on the true mass hierarchy. 
The discovery reach of the experiment is therefore expected to be dependent on the true mass 
hierarchy. Likewise, while the value of 5cp in the fit is inconsequential as sin^ 26*13 = in the fit, 
the data itself and hence the discovery reach, would depend on (5cp(true). For each sin^ 2^13 (true), 
we generate the data for all possible values of 5cp(true) and for both the mass hierarchies. For each 
case, the data is then fitted assuming sin^ 2^13 = and marginalizing over the other oscillation 
parameters, returning a value of Xmin f*^^ each data set. We choose the minimum amongst these 
Xmin fi^*^ value of sin^ 2^13 (true) for which we could claim a signal in the detector at the 
3(7 C.L. In the right panel of Fig. 10 we show this "most conservative"^^ sin^ 2^i3 discovery reach 
of our experiment as a function of 7. We assume equal 7 for both the ions for the red solid curve. 
One can see that for 7 = 650, the most conservative discovery reach is sin^ 2^i3(true) = 5.11 x 10"'^ 
while if (5cp(true) = then this will be sin^ 2^13 (true) = 5.05 x 10""^. For the blue dashed line 
we assume that the 7 for ^Li is scaled down by a factor of 1.67 compared to that for ^B, plotted 
on the X-axis. Since for same 7, neutrino is the statistically weaker channel, the red line mainly 
corresponds to what we expect for the true NH. For the scaled 7 case since the antineutrino 
channel becomes statistically weaker, the lower comes from this channel and the blue dashed 
line corresponds to what we expect for the true IH. 

4.3 sin^ 26*13 Precision 

In Fig. 11 we show how precisely the mixing angle sin^2^i3 will be measured, if we observe 
a 6^13 driven signal at the detector. The left panel depicts as a function of sin^ 26^13 (true) the 
corresponding range of allowed values of sin^ 26*13 at the Scr C.L. We have assumed 7 = 500 and 
5cp — The sohd line is assuming NH to be true, while the dashed fine is for IH true. Note that 
in the fit we always marginalize over the hierarchy and 5cp- The right panel shows the variable 

^^This is "most conservative" in the sense that no matter what the choices of (5cp(true) and the true neutrino 
mass ordering, the ^13 discovery hmit cannot be worse than the value presented. 
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Figure 11: The precision with which sin^ 26'i3 will be measured by the CERN-INO beta-beam 
experiment as a function of sin^ 26'i3(true). Left panel shows the 3cr allowed range of sin^2^i3 
while the right panel shows the precision defined in the text. 



"precision" which we define as 



(sin 



(sin^ 2^,3)_ + (sin^ 2^,3)^ ' 

where (sin^ '^(^n)max (^^^^ '^^^3)min maximum and minimum allowed values of sin^ 26*13 

respectively at 3a. 



5 Conclusions 

Long baseline experiments which will use the golden P^^ channel for determining the neutrino 
oscillation parameters face a serious threat from the menace of clone solutions due to the so- 
called parameter degeneracies. These degeneracies come in three forms: the Sep — 613 intrinsic 
degeneracy, the 6cp — sgn{Aml^) degeneracy and the 6*23 octant degeneracy, and necessarily result 
in degrading the sensitivity of the experiment. The CERN-INO near-magic distance of 7152 km 
offers the possibility of setting up an experiment at a baseline where the 6c p dependent terms 
almost drop out from the expression of the golden channel probability. Thus two out of the 
three degeneracies are evaded, providing a platform for clean measurement of 6*13 and sgn{Amli), 
two major players in our understanding of the origin of neutrino masses and mixing. A large 
magnetized iron calorimeter with a total mass of at least 50 kton is expected to be built soon at 
INO. It will be ideal for detecting multi-GeV z/^ and hence can be used as the far detector for a 
high energy beta-beam. 

In this paper we studied in detail the physics reach of the CERN-INO magical beta-beam 
set-up and extended our analysis presented in [19]. Most importantly, we showed the impact of 
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adding data from both the neutrino and antineutrino runs of the experiment. Combined data 
from both the neutrino and antineutrino polarities brings a major improvement in the hierarchy 
sensitivity of the experiment. We took into account the uncertainty due the solar parameters 
Am|]^ and sin^ 612 and marginalized the over them. We also accounted for the uncertainty 
in the Earth matter density profile in our fit- We studied the impact of changing the energy 
threshold and the background rejection factor of the detector. We probed the importance of using 
the full spectral information on the final sensitivity of the experiment. For 7 = 650, 5c;p(true) = 
and true NH, the sensitivity to hierarchy determination at 3a improves almost two orders of 
magnitude from sin^ 2^13 (true) = 1.15 x 10~^ for the neutrino channel using the total rate only to 
sin^ 2^^i3(true) = 5.51 x 10~^ when full spectral data from neutrino and antineutrino channels are 
combined. Even though the effect of 5cp (true) on the event rate of our experiment is expected 
to be small, there is some residual dependence on it because the CERN-INO distance does not 
conform to the exact magic baseline. We studied the change in the hierarchy sensitivity due to the 
uncertainty in Sep- It turns out that for 7 = 650 and with NH (IH) true, the best sensitivity to 
hierarchy determination corresponds to sin^ 2^13 (true) = 3.96x10^"' (sin^ 2^^i3(true) = 2.96x10"*^), 
while the worst case is sin^ 2^i3(true) = 5.58 x 10"^ (sin^ 2^i3(true) = 3.59 x 10"^). 

We presented a detailed analysis of the potential of probing 6*13 at this experiment. We defined 
and studied the ^13 reach in terms of three performance indicators: the sensitivity reach, the 
discovery reach and the precision of sin^ 2^13 measurement. The sensitivity reach is defined as the 
upper limit on sin^ 2^13 that we would be able to impose in case the data is consistent with null 
oscillations. At 3a C.L. the sensitivity reach corresponds to sin^ 26*13 = 1.14 x 10~^ for 7 = 650 and 
this is independent of the true hierarchy and 5cp (true). The discovery reach is defined as the true 
value of the mixing angle for which we have an unambiguous oscillation signal in the detector. 
At 3a C.L. the discovery reach corresponds to sin^ 2^13 (true) = 5.05 x 10~^ (sin^ 2^13 (true) = 
2.96 X lO^'') for 7 = 650, 5cp(truc) = and NH (IH) true while the most conservative limit 
irrespective of 5c;p(true) and the true neutrino mass ordering is sin^ 26'i3(true) = 5.11 x 10"*^. We 
also presented the expected precision with which sin^ 2^13 would be determined in this experiment 
for sin^ 2^i3(true) > 10"^ 

Neutrino physics is in wait for the next great leap forward in the decade ahead. Beta-beams and 
an iron calorimeter detector at a very long baseline may well turn out, as we have demonstrated, 
to be a key player in this endeavour. 
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Appendix: The Numerical Analysis 

In this appendix, we give more details of our numerical technique through which we statistically 
explore the physics potential of our set up. The induced ii~ event spectrum at INO is estimated 
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using 



N,^T Tin fine / dE dEA(t>{E) a,^{E) R{E,Ea) Pe^{E) (5) 

where are the number of events in energy bin i with lower energy hmit and upper energy 
Umit E™^^^ T is the exposure time, n„ are the number of target nucleons, is the charge 
identification efficiency, e is the detection efficiency, (j){E) is the beta-beam flux at INO, a^j^ is the 
detection cross section for i/^, R{E,Ea) is the detector energy resolution function^^, E being the 
true energy of the incoming neutrino and Ea the measured energy of the muon. The corresponding 

spectrum due to is given by replacing by a^,^ by ap^ and Pg/^ by Pgp, in Eq. (5). 
We have used fjo — 0.95, e = 0.8 and taken the neutrino-nucleon interaction cross sections from 
[46, 47, 48]. In this paper we consider only the neutral current backgrounds coming from the beam 
itself. We assume that the background can be rejected very efficiently by imposing suitable cuts 
such that only a fraction of 10~^ of these neutral current backgrounds survive. The background 
is assumed to have the same shape as the signal. But one should keep in mind the fact that this 
shape is not much of an issue since anyway the background is very small. 

For our statistical analysis we employ a function defined as 

2 2 _|_ 2 _|_ 2 //^\ 

where the first term is the contribution from the neutrino channel, the second term comes from the 
antineutrino channel, while the last term comes from imposing priors on the oscillation parameters 
which we allow to vary freely in our fit and which we expect will be determined better from other 
experiments at the time when the data from the CERN-INO beta-beam set-up would be finally 
available. The fo^ f^e neutrino channel is given by 



max 



where n is the total number of bins. 



I i=l 



(7) 



UM, {^s, 6}) = ATf ({u;}) [1 + TT^e,] + N^[1 + 7r%\ , (8) 

Nj^{{u!}) given by Eq. (5) being the predicted number of events in the energy bin i for a set of 
oscillation parameters cu and are the number of background events in bin i. The quantities tt* 
and TT^ in Eq. (8) are the systematical errors on signals and backgrounds respectively. We have 
taken tt* = 2.5% and tt^ = 5%. The quantities and are the "pulls" due to the systematical 
error on signal and background respectively. The data in Eq. (7) enters through the variables 
Xi = N^^ + N^, where Nf^ are the number of observed signal events in the detector and iVf is the 
background, as mentioned earlier. We simulate the signal event spectrum using Eq. (5) for our 
assumed true values for the set of oscillation parameters. These assumed true values are given in 
the first column of Table 2. For sin^ 2^13 (true), Sep (true) and the true hierarchy we use different 



^^We assume a Gaussian resolution function with a = 0.15E. 
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|Am2^(true)| = 2.5 x lO'^ eV^ 




1.5% 


8111^2^23 (true) = 1.0 


(7(sin2 2^23) 


= 1% 


Amii(true) = 8.0 x 10"^ eV^ 


(7(Amii) = 


2% 


sin2^i2(true) = 0.31 


(T(sin^ 012) = 


= 6% 


p(true) = 1 (PREM) 


a{p) = 5% 





Table 2: Chosen benchmark values of oscillation parameters and their la estimated error. The 
last row gives the corresponding values for the Earth matter density. 



options and mention them wherever applicable. In our fit we marginalize over all oscillation 
parameters, the Earth matter density, as well as the neutrino mass hierarchy, as applicable. We 
do this by allowing all of these to vary freely in the fit and picking the smallest value for the 
function. Of course we expect better determination of some of these parameters, which are poorly 
constrained by this experimental set-up. Therefore, we impose a "prior" on these parameters 
through the xLior given by 



2 ^1 J 1 ^ I , 

A prior I _ / 1 A 2 I \ I 



I Amgi I - I Amgi (true) | V / sin^ 26*23 - sin^ 26*23 (true) " 



<^prior 



+ 



(7(|Amii|) J \ (7(sin2 2e23) 

/ Am|i - Am|i (true) V / sin^ 6*12 - sin^ 6*12 (true) " 

(7(Amli) J ^ [ a{sm^ e^) . 

2 



- m ■ 

where the la error on these that we use are taken from [16, 49] and are given in the left column 
of Table 2. In our computation, we have used a matter profile inside the Earth with 24 layers. In 
Eq. (9), p is a constant number by which the matter density of each layer has been scaled. The 
external information on p is assumed to come from the study of the tomography of the earth [50] . 
In Eq. (9), p varies from 0.95 to 1.05 i.e., 5% fiuctuation around 1. 

Note that in our definition of the function given by Eq. (6) and (7) we have assumed that 
the neutrino and antineutrino channel are completely uncorrelated, all the energy bins for a given 
channel are fully correlated, and and are fully uncorrelated. We minimize the Xtotai two 
stages. First it is minimized with respect to and ^5 to get Eq. (7), and then with respect to 
the oscillation parameters uj to get the global best-fit. For minima with respect to and we 
require that 

^ = and ^ = . (10) 
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From Eq. (7, 8, 10) we get, 



On ai2 

^21 0,22 



6 



Cl 
C2 



where. 



C2 = E 



i=l 
n 

022 = 



i=l 



{Nl^ + ATf )2 



+ 1 , 



(iVf + iVi')2 

n 

Ol2 — O21 — Y, 



+ 1 , 



i=l 



(11) 



(12) 



Using Eq. (11 and 12), we calculate the values of and ^5 and then we use these values to 
calculate xl^^u ■ 1^ ^ similar fashion, we estimate xl^^p to obtain the Xtotai- 
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